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Abstract 

We show that in supersymmetric extensions of the Standard Model gluino box 
diagrams can yield a large A/ = | contribution to s — > dqq FCNC processes, which 
may induce a sizable CP-violating contribution to the 1 = 2 isospin amplitude in 
K — > vrvr decays. This contribution only requires moderate mass splitting between 
the right-handed squarks ur and da, and persists for squark masses of order 1 TeV. 
Taking into account current bounds on Imd^J" from K-K mixing, the resulting 
contribution to e'/e could be an order of magnitude larger than the measured value. 
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The recent confirmation of direct CP violation in K tttt decays is an impor- 
tant step in testing tlie Cabibbo-Kobayaslii-Maskawa (CKM) meclianism for CP vio- 
lation in the Standard Model. Combining the recent measurements by the KTeV and 
NA48 experiments with earlier results from NA31 and E731 gives Re(e'/e) = 
(2.12 ± 0.46) X 10~^. This value tends to be higher than theoretical predictions in the 
Standard Model, which center below or around 1 x 10^'^ 0. Unfortunately, it is difficult 
to gauge the accuracy of these predictions, because they depend on hadronic matrix 
elements which at present cannot be computed from first principles. A Standard-Model 
explanation of e'/e can therefore not be excluded. Nevertheless, it is interesting to ask 
how large e'/e could be in extensions of the Standard Model. 

In the context of supersymmetric models, it has been known for some time that it 
is possible to obtain a large contribution to e'/e via the AI = ^ chromomagnetic dipole 
operator without violating constraints from K-K mixing . It has recently been pointed 
out that this mechanism can naturally be realized in various supersymmetric scenarios 
0. In this Letter we propose a new mechanism involving a supersymmetric contribution 
to e'/e induced by A/ = | penguin operators. These operators are potentially important 
because their effect is enhanced by the A/ = ^ selection rule. Unlike previous proposals, 
which involve left-right down-squark mass insertions, our effect relies on the left-left 
insertion 5^^ and requires (moderate) isospin violation in the right-handed squark sector. 

The ratio e'/e parametrizing the strength of direct CP violation in i^T — >■ tttt decays 
can be expressed as 

e y2|e| VReAs ReAo) ' ^ ' 

where Aj are the isospin amplitudes for the decays —>■ (vrvr)/, 5/ are the corresponding 
strong-interaction phases, and the ratio u = IieA2/IieAo ^ 0.045 signals the strong 
enhancement of A J = | transitions over those with AI = |. From experiment, we take 
|e| = (2.280 ± 0.013) x 10^^ and 0, = (43.49 ± 0.08)° for the magnitude and phase of 
the parameter e measuring CP violation in K-K mixing, and ^2 — 5o = —(42 ±4)° for 
the difference of the S-wave vrvr scattering phases in the two isospin channels. It follows 
that, to an excellent approximation, e'/e is real. 

In the Standard Model, the isospin amplitudes Aj receive small, CP-violating con- 
tributions via the ratio iYts^td) / {y*gVud) of CKM matrix elements. This ratio enters 
through the interference of the s uud tree diagram with penguin diagrams involving 
the exchange of a virtual top quark. According to (P, contributions to e'/e due to the 
AI = I amplitude Im742 are enhanced relative to those due to the AI = | amplitude 
ImAo by a factor of ~ 22. However, in the Standard Model the dominant CP- 
violating contributions to e'/e are due to QCD penguin operators, which only contribute 
to Aq. Penguin contributions to A2 arise through electroweak interactions and are sup- 
pressed by a power of a. Their effects on e'/e are subleading and of the same order as 
isospin-violating effects such as 'K^-rj-rj' mixing. 

Here we point out that in supersymmetric extensions of the Standard Model poten- 
tially large, CP-violating contributions can arise from fiavor-changing strong-interaction 
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processes induced by gluino box diagrams. Whereas in the hmit of exact isospin symme- 
try in the squark sector these graphs only induce AI = | operators at low energies, in 
the presence of even moderate up-down squark mass splitting they can lead to operators 
with large AI = | components. In the terminology of the standard effective weak Hamil- 
tonian, this implies that the supersymmetric contributions to the Wilson coefficients of 
QCD and electroweak penguin operators can be of the same order. Specifically, both sets 
of coefficients scale like a^/m^ with m a generic supersymmetric mass, compared with 
agaw/m"^ and aavy/m^, respectively, in the Standard Model. These contributions can 
be much larger than the electroweak penguins of the Standard Model even for super- 
symmetric masses of order 1 TeV. On the other hand, supersymmetric contributions to 
the Wilson coefficients proportional to electroweak gauge couplings are parametrically 
suppressed and will not be considered here. 

We find that the relevant AS" = 1 gluino box diagrams lead to the effective Hamilto- 
nian 

G ^ ~ 
KflF = ^ ^ k(^) + cf (/i) (/i)] + h.c. , 

V ^ i=l 

where 

Q\ = idaSa)v-A iQl3q(i)v+A , 

Q2 = idaSi3)v~A {Qf3<la)v+A , 

Ql = idaSa)v-A {.(Il3<ll5)v-A , 

Q\ = idaSfj)v-A {.(Il3<la)v-A 

are local four-quark operators renormalized at a scale fi <^rn, Qf are operators of oppo- 
site chirally obtained by interchanging V—A ^ V+A, and a summation over q = u,d, . . . 
and over color indices a, j3 is implied. In the calculation of the coefficient functions we 
use the mass insertion approximation, in which case the gluino-quark-squark couplings 
are flavor diagonal. Flavor mixing is due to small deviations from squark-mass degen- 
eracy and is parametrized by dimensionless quantities 5^^, where z,j are squark flavor 
indices and A,B refer to the chiralities of the corresponding quarks (see, e.g., In 
general, these mass insertions can carry new CP-violating phases. Contributions in- 
volving left-right squark mixing are neglected, since they are quadratic in small mass 
insertion parameters, i.e., proportional to S^^dq^. We define the dimensionless ratios 

\ m-g ) \ m-g J 

where mnj^ ^ and "^j^ ^ denote the average left- or right-handed squark masses in the up 
and down sector, respectively. S'U{2)l gauge symmetry implies that the mass splitting 
between the left-handed up- and down-squarks must be tiny; however, we will not assume 
such a degeneracy between the right-handed squarks. For the Wilson coefficients cj at 
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the super symmetric matching scale m we then obtain 



where 
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The results for the coefficients c- are obtained by interchanging L <-> in the expressions 
for cj. 

It is straightforward to relate the quantities cf to the Wilson coefficients appearing 
in the effective weak Hamiltonian of the Standard Model as defined, e.g., in 0. We find 
(At = V*Vtd) 



C3 + 2cf 
3 

+ 24 



(-At) C, = ^1±^ 

(-At) a = '-^±^ 



for the QCD penguin coefficients, and 

3 

-(— Ai) = c 



C? = Ac,; 



1...4 



(2) 



for the coefficients of the electroweak penguin operators. The supersymmetric contri- 
butions to the electroweak penguin coefficients vanish in the limit of universal squark 
masses. However, for moderate up-down squark mass splitting the differences Aq = 
c" — cf are of the same order as the coefficients cf themselves. In this case gluino 
box contributions to QCD and electroweak penguin operators are of similar magnitude. 
This conclusion is unaltered when additional contributions to C3...6 from QCD penguin 
diagrams with gluino loops are taken into account. Because the electroweak penguin 
operators contain A/ = | components their contributions to e'/ e are strongly enhanced 
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and thus are expected to be an order of magnitude larger than the contributions from the 
QCD penguin operators. In this Letter, we focus only on these enhanced contributions. 

The renormalization-group evolution of the coefficients Aq (and Acj) from the super- 
symmetric matching scale m down to low energies is well known. In leading logarithmic 
approximation, one obtains B 



Ac2(m) 



Ac,(;z) + ^ = .«/^o 



Aci(m) 



Acsifi) + Ac4(/i) = [Ac3(m) + Ac4(m)] , 

Acsifi) - Ac4(/i) = K^/^" [Ac3(m) - Ac4(m)] , 

where k = as(/i)/as(m), and Po = 11 — ^nj. It is understood that the value of Pq is 
changed at each quark threshold. We use the two-loop running coupling normalized to 
«s(^z) = 0.119 and take the quark thresholds at rrit = 165 GeV, = 4.25 GeV and 
rric = 1.3 GeV. In Table ^ we give the imaginary parts of the coefficients Aci^2 and Ac3^4 at 
the scale /i = rric, obtained for the illustrative choice m = rrig = m^^ = mj^ = 500 GeV 
and three (larger) values of m^^. Since the mass splitting between the left-handed ul 
and (11 squarks is tiny, we can safely neglect the coefficients Ac3,4 and Aci,2 in our 
analysis. Note that for fixed ratios of the supersymmetric masses the values of the 
coefficients scale like m~^, i.e., significantly larger values could be obtained for smaller 
masses. For comparison, the last column contains the imaginary parts of Aci...4 in the 
Standard Model computed from (0) using ImA^ = 1.2 x 10""^ and the next-to-leading 
order Wilson coefficients Cj compiled in . We observe that for supersymmetric masses 
of order 500 GeV, and for a mass insertion parameter lm6^j" of order a few times 10~^ 
(see below), the Wilson coefficient Ac2 can be significantly larger than the value of the 
corresponding coefficient in the Standard Model, which is proportional to Cg. This is 
interesting, since even in the Standard Model the contribution of Cg to e'/e is significant. 

In estimating the supersymmetric contribution to e'/e we focus only on the {V — A)^ 
{V + A) operators associated with the coefficients Aci and Ac2, because their matrix 
elements are chirally enhanced with respect to those of the {V — A)®{V — A) operators. 
The penguin operators contribute to the imaginary part of the isospin amplitude A2. 
The real part is, to an excellent approximation, given by the matrix elements of the 
standard current-current operators in the effective weak Hamiltonian. Evaluating the 
matrix elements of the four-quark operators in the factorization approximation, and 
parametrizing nonfactorizable corrections by hadronic parameters as defined in 
we obtain 

ImAr'' 3 ml lm[Ac2{m,) + \Aci{mc)] (m^) 



ReA2 2 m2(m,) - ml{m^) z+{m,) B?\ 



rrir 



Following common practice we have neglected a tiny contribution proportional to the 

(2) (2) 

difference Bj — B^ . In the above formula 2+ is a combination of Wilson coefficients. 
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Table 1: Values of the coefficients Aci{mc) and Acj(mc) in units of 10~^ Im S^J" 
and 10~^Im5^^, respectively, for common gluino and down-squark masses of 
500 GeV and different values of rriu^. The last column shows the corresponding 
values in the Standard Model in units of 10~^. 
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>c) 


0.56 


0.87 


1.17 


-7.63 



The product z+ = 0.363 is scheme independent and can be extracted from experi- 
ment. Note that at leading logarithmic order the scale dependence of the combination 
Ac2 + |Aci cancels the scale dependence of the running quark masses, and hence the 

(2) 

hadronic parameter is scale independent. 

Putting everything together, we find for the supersymmetric A/ 



I contribution to 



19.2 
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where 



32 2 

X(x, y,z) = — [/(x, y) - /(x, z)] 



27 ^ ^ V ^ /J 27 
The existence of this contribution requires a new CP-violating phase 0^ defined by 
= 1 5^/1 sin 02.- The measured values of Aitlk and e in K-K mixing give bounds 
on Re and Im(5^/')^, respectively, which can be combined to obtain a bound on 

Im^f/' as a function of (pi- Using the most recent analysis of supersymmetric contribu- 
tions to K-K mixing in , we show in Figure |I| the results obtained for m^^ = 500 GeV 
and three choices of m^. 

It is evident that the bound on \m.5^j^ depends strongly on the precise value of 
(pL. To address the issue of how large a supersymmetric contribution to e'/e one can 
reasonably expect via the mechanism proposed in this Letter, it appears unnatural to 
take the absolute maximum of the bound given the peaked nature of the curves. To 
be conservative we evaluate our result (0) taking for Im 5^^ one quarter of the maximal 
allowed values shown in Figure |l], noting however that a larger effect could be obtained 
for the special case of \(f)i\ very close to 90°. Our results for |e'/e| are shown in Figure ^ 
as a function of m^^ for the case mj^ = m^^ = 500 GeV and the same three values 
of Trig considered in the previous figure. The choice mj^ = mj^ is made for simplicity 
only and does not affect our conclusions in a qualitative way. Except for the special 
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Figure 1: Upper bound on |Im(5f/| versus the weak phase (in degrees) for 
= SOOGeVand {jng/m^^Y = 1 (sohd), 0.3 (dashed) and 4 (short-dashed). 

case of highly degenerate right-handed up- and down-squark masses, the AI = | gluino 
box-diagram contribution to e'/e can by far exceed the experimental result, even taking 
into account the bounds from Atjik and e. Indeed, even for moderate splitting Figure || 
implies substantially stronger bounds on |Im5f/'| than those obtained from K-K mixing. 
This finding is in contrast to the commonly held view that supersymmetric contributions 
to the electroweak penguin operators have a negligible impact on e'/e. In this context, it 
is worth noting that a large mass splitting between and can be obtained, e.g., in 
GUT theories without SU(2)j:j symmetry and with hypercharge embedded in the unified 
gauge group, without encountering difficulties with naturalness 

The allowed contribution to e'/ e in Figure |^ increases with the gluino mass (for fixed 
squark masses) because the K-K bounds become weaker in this case. If all supersym- 
metric masses are rescaled by a common factor ^, and the bounds from K-K mixing are 
rescaled accordingly, the values for e'/e scale like 1/^ modulo logarithmic effects from the 
running coupling ^^(^m). Therefore, even for larger squark masses of order 1 TeV the 
new contribution to e'/e can exceed the experimental value by a large amount, implying 
nontrivial constraints on Im^f/". 

Before concluding, we note that in the above discussion we have made no assumption 
regarding the mass insertion parameter lm.5fj^ = sin0^ for right-handed squarks. 
In models where \S^J^\ is not highly suppressed relative to l^f/" |, much tighter constraints 
on Im 6^J" can be derived by applying the severe bounds on the product 6^J^ obtained 
from the chirally-enhanced contributions to K-K mixing. In analogy with Figure 0, we 
obtain an upper bound on llm^f/^l as a function of (pi and which is sharply peaked 
along the line (j)L + (pR = mod tt and scales like \5gJ" /5^^\^^'^. As above, we take one 
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Figure 2: Supersymmetric contribution to |e'/e| (in units of 10~^) versus 
rriuj^, for ms{mc) = 130 MeV, Bg {rric) = 1, m^^ = m^^ = 500 GeV, and 
{rrig/m^^Y = 1 (solid), 0.3 (dashed) and 4 (short-dashed). The values of 
|Im5^/'| corresponding to the three curves are 0.011, 0.005 and 0.027, respec- 
tively (see text). The band shows the average experimental value. 

quarter of the peak value obtained using the results compiled in [0 . Considering the case 
mj^ = nig = 500 GeV for example, we find that the upper bounds on |Im5f;f | a reduced 
by a factor ranging from 3% in the limit where \SfJ^\ = to 8% for \SfJ^\ = 0.1 1 5^/' |. 
In the latter case, the supersymmetric contribution to e'/e can still be of order 10~^, i.e., 
comparable to the measured value. Moreover, significantly larger values can be obtained 
for special points in moduli space, where the weak phases obey 0l + 0i? ~ mod vr. 

In summary, we have shown that in supersymmetric extensions of the Standard Model 
gluino box diagrams can yield a large AI = | contribution to e'/e, which only requires 
moderate mass splitting between the right-handed squarks, i.e., {'rriuB —^dj^/^dji ^ 
In a large region of parameter space, the measured value of e'/e implies a significantly 
stronger bound on Im 5^^ than is obtained from K-K mixing. 
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